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Abstract: Recent data for the ratio R{Q) = QF2{Q'^)/ Fi{Q'^) shocked the 
community by disobeying expectations held for 50 years. We examine the sta- 
tus of perturbative QCD predictions for hehcity-flip form factors. Contrary 
to common behef, we find there is no rule of hadron helicity conservation for 
form factors. Instead the analysis yields an inequality that the leading power 
of helicity-flip processes may equal or exceed the power of helicity conserving 
processes. Numerical calculations support the rule, and extend the result 
to the regime of laboratory momentum transfer Q^. Quark orbital angular 
momentum, an important feature of the helicity flip processes, may play a 
role in all form factors at large Q^, depending on the quark wave functions. 



1 The Nature of High Energy Reactions 



There is an asymmetry in high energy reactions due to the Lorentz trans- 
formation. The spatial coordinate parallel to the boost axis is Lorentz con- 
tracted. The momentum fraction x of partons inside hadrons is thereby 
distributed over the entire possible range < x < 1. This phenomenon 
is dynamical, because a boost of interacting fields is dynamical, and the x 
dependence of wave functions cannot reliably be calculated in perturbative 
QCD. Instead, the x-dependence of wave functions is extracted from experi- 
mental data. 

Meanwhile the transverse spatial coordinate h is Lorentz invariant. The 
transverse coordinate has a certain calculability via perturbative QCD. There 
is a great deal of interest and controversy associated with the transverse co- 
ordinate in exclusive reactions. The transverse spatial coordinate can be 
probed in reactions sensitive to the angular momentum flow. In some reac- 
tions, the sum of the helicities going into a hadronic reaction is automati- 
cally conserved. This is the case of the proton's Dirac form factor Fi in the 
high energy limit. When the sum of the helicities is not conserved, angular 
momentum conservation requires either extra constituents, or quark orbital 
angular momentum. This is the case of F2, the proton's Pauli form factor. 

The Jefferson Laboratory has observed Q F2{Q'^) / Fi{Q^) ~ const, up to 
the highest values of ~ 5.8 GeV^ yet measured [H |2]. The data was 
initially very surprising, and the field may have reached a pivotal point in 
comparison with the quark model. For a long time it was held sacred that 
Q'^F2/Fi ~ const at large Q^. This rule appears to predate QCD. It has 
ancient origins in renormalization questions involving protons as elementary 
fundamental fields jS]- 

Meanwhile a perturbative QCD model assuming non-zero quark orbital 
angular momentum (0AM) predicted Q F2{Q^) / Fi{Q^) ~ const ||1 E]. A 
relativistic quark model prediction fits the flatness of Q F2{Q^) / Fi{Q'^) 
equally well. These papers countered the ancient wisdom, because they 
shared the common feature of quark orbital angular momentum in the wave 
functions. 

The hypothesis of zero 0AM sometimes appeals to the non-relativistic 
quark model in the rest frame. Yet very little of high energy physics and 
pQCD starts in the rest frame. The perturbative quark wave functions are 
unrestricted in angular momentum content, except for Lorentz symmetry. 
Observation of non-zero quark 0AM is a leading candidate to resolve the 
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proton "spin crisis" , which is the fact that the sum of the perturbative quark 
spins does not equal the spin of the proton. Consequently the Jefferson lab's 
data has even broader implications than the mystery of large form factors. 

Here we address whether Q F2{Q'^) / FiiQ"^) ~ const, is a transient feature 
of comparatively low energy experiments, or a fact destined to persist at 
higher Q^. If the flatness of this ratio is due to quark 0AM, will the ratio 
stay fiat with increasing Q^l The question leads us to re-examine the roots of 
the "hadron helicity conservation" rule [7j. We find that pQCD itself is rather 
neutral, and only gives an inequality between the powers governing helicity- 
fiip processes and helicity conserving ones. As a result QF2/ Fi ~ const, may 
extend to arbitrarily large values of Q^, without violating anything sacred. 

2 Definitions and Discussion 

Orbital angular momentum: By quark orbital angular momentum (0AM) 
we refer to an expansion in S0{2) representations (commonly known as 
states) with quantization axis aligned along the particle 3-momentum P. 
The eigenstates of are invariant under boosts along the z-axis. We do not 
use spherical harmonics, and we treat the longitudinal coordinates as scaling 
variables. Let the transverse spatial coordinate be 6 = 6 (cos0, sincj)). We 
expand operators or wave functions as 

V' = Ee'"Vm(x, b) 

m 

where x is the Feynman light-cone fraction of z (or "+") momentum.^ 0AM 
can also be re-expressed with generalized parton distributions in a manifestly 
gauge-invariant formalism |Hl E] . 

Chirality versus Helicity: The proton's Pauli form factor F2 contains 
information on the orbital angular momentum of the quarks, but it is indirect. 
Strictly speaking, the amplitude iu{p', s')F2{Q'^)a^^qi,u{p, s) represents the 
amplitude for chirality of the proton to flip under momentum transfer Q. 
The chirality (eigenvalue of 75) of light quarks flips very little in pQCD 

^Of course perturbative wave functions are theory constructs, and fields at finite sepa- 
ration are technically not gauge invariant, in general. Observables are nevertheless gauge 
invariant when all the legs of diagrams are contracted properly. 
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reactions, by about nig/Q « 1, for quarks of mass niq ~ fewMeV. The 
chirality of Fermions is proportional to their hehcity, with corrections of order 
m/Q. Adding another constituent to the scattering is also down by 0(1/Q^). 
Then by this reasoning, at large > GeV^, it is not possible to flip either 
the chirality or helicity of the proton with a virtual photon, unless there 
is internal quark orbital angular momentum present to satisfy the selection 
rules. 

What is Quark Counting? We separate the quark-counting model of 
Brodsky, Farrar, and independently Matveev et al ^0] from the asymptotic 
short distance [asd) model of Brodsky and Lepage fUE]. The earlier theory 
is one of counting propagators and the number of scattered constituents. 
The latter theory is much more detailed, imposing a certain factorization 
of the hard reaction into components made from the s-wave, m = Bethe- 
Salpeter wave functions. This does not come by listing all diagrams, and 
is not a feature of the starting theory. Instead the framework is developed 
by assuming the framework and classifying terms within the assumptions. 
The asd approach is characterized by taking the zero-distance limit in the 
first step, and replacing the rest of the calculation by integrals over Feynman 
X fractions using "distribution amplitudes", or similar quantities with no 
transverse information. However the Feynman rules instruct us to leave the 
longitudinal and transverse integrals coupled. Perform the integrals, and 
afterwards take the limit of large Q (if wanted). If the two limits are not the 
same (and they are not in general), then the asd assumptions can fail. 

The transverse coordinates are gone in the asd approach, being evalu- 
ated within 1/Q of zero. HHC follows instantly as a test of the framework, 
independent of the wave functions used. By omitting 0AM m ^ 0, the pre- 
diction for F2 = 0, and one cannot recover any prediction for finite F2 in 
the formalism. As far as we know, all asd predictions of F2 are indirect and 
deduced by elimination: since Fi ~ l/Q"^, and QF2 could not be calculated, 
then F2 ~ Fi/Q'^ must lie in the detritus not calculated. A recent calcu- 
lation however, finds that within the framework of a generalized asd 
model F2 ~ l/<5^- It is also possible to modify the model by including quark 
mass effects [21 • Our goal here is to understand the limit of arbitrarily light 
quarks and neglecting effects of order rUq/Q. 

We separated quark counting from asd because the two are not the same 
theory. Should one believe either? There are good indications from the 
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scaling laws that quark-counting has some truth. We have no religious com- 
mitment here, and the possibility that a fraction of the amplitude is due 
to "soft physics" must be given credit [13 Cn] ■ But amplitudes cannot all 
be soft, because the form factors are not seen to fall exponentially with Q. 
Indeed the quark-counting scaling laws generally work well. Meanwhile the 
inapplicability of HHC to F2 or any other helicity flip reactions ^7j shows 
we cannot use asd. 

Regarding Generalized Parton Distributions: Generalized parton dis- 
tributions (GPD) appear to be an ideal way to proceed. Lorentz covariance 
can be used to set up matrix elements and expressions for the form factor in 
an apparently model-independent way. As far as we know, form factors were 
the first instance of GPD, used in the paper of Soper ^H] in 1977, which also 
contains a transformation to a particular transverse spatial coordinate. We 
will have occasion to revisit the conclusions of that paper in Section 3.7. 

Despite our GPD-hased predictions jSl IH EI for ratios and the welcome 
rediscovery of GPD in the field, we chose not to make them the vehicle for 
this analysis. The reason is that GPD are so general they do not immediately 
contain the information there are three quarks in the proton. To incorporate 
the information one can start with wave functions and integrate out all but 
two quark legs. Since our concern is precisely these integrations, we would 
have nothing to gain. We caution the reader, in any event, that the method 
of approximating integrals by asd methods will have similar limit- interchange 
problems when clothed in GPD language. Diehl et al. [19 discuss general 
relations, and their formulas codify a relation of F2 to quark 0AM. 

Some Familiar, But Inexact Assertions: Suppose we have non-zero 
0AM in the wave functions and we try to use the assertion of the asd ap- 
proach, 

A6AQt > 1 • (1) 

Here A6 is the resolved transverse quark separation, in a frame where the 
momentum transfer Q ~ AQt is transverse. We have written the relation like 
the uncertainty principle, to give it a chance to be seductive. We observe next 
that, merely from continuity, a wave function ipm carrying m units of angular 
momentum scales like 6*" as 6 — > 0. Then under these assumptions each unit 
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of orbital angular momentum of the quarks will lead to amplitudes suppressed 
by a corresponding power of 1/Q at large Q. This familar assertion has been 
repeated endlessly in the literature, yet we will show that this type of counting 
does not represent QCD. It is seductive but it is not right. 

Our Approach: To address F2 properly, one must go beyond short-distance 
to restore the transverse coordinate. This is called "impact-parameter factor- 
ization" I21j. This well-justified method has dominated recent attention 
in perturbative QCD [221 123 1211 12SI ISHl • In impact-parameter factorization 
the form factor is written as 

F(g2) = J{Udx,d^kT,i) {Udx'id^k'T^i) ^{x,kT)H{xi,x'i,kT,i,k!r/, Q)^(x', A;^). (2) 

The impulse approximation has been used to set the light-cone "time" to 
zero. There are no a-priori assumptions in Eq. |21 about short-distance. If 
Sudakov effects are used consistently, then wave functions in Eq. |21 concen- 
trate the dominant region into one which is perturbatively calculable, without 
assuming zero- distance as a starting point. We do not go to the further ex- 
treme of Ref. [2^ towards asserting that the ultimate output of Sudakov 
effects is the asd model. We find this unjustified. To get that result one 
needs assumptions about the wave functions that are simply not known. 

3 Calculations 

Here we present pQCD calculations and also illustrate certain features that 
are even more general. 

3.1 Role of The Transverse Coordinate 

Return to the general expression Eq. |21 To leading order and neglecting 
transverse momentum in Fermion numerators, the hard scattering depends 
on differences fcy — k'j, before and after the hard collision. The importance 
of this variable seems rather general, because the sum of the transverse mo- 
menta are conjugate to the overall spatial location of the hard scattering, 
which by translational invariance drops out.^ In a process with a single hard 

^From translational invariance, the hard kernel must depend on the difference of space 
coordinates. The transverse separation b is only one such difference, and the one of interest 
here. Other differences such as the longitudinal ones also occur. 
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exchange (a pion or meson form factor) the transverse integrals take the form 
of convolution: 

F{Q) = J d'kTfikTmkr); 
f{kT) = J Sk'^H^kr - k'^] QM^'t), 

The X variables were suppressed. Consequently the expressions are diagonal 
in b, with 

F{Q) = I ^^dxdx'rib, x)H{b;x, x', Q)tPib, x'), 
= E I -^dxdx'dcPe-'^"'-"'''^'f'r^ib,x)H„,ib;^^ 

All contributions to 0AM are explicit at this stage. Expansion of the hard 
scattering was introduced because the kernel can also carry angular mo- 
mentum and be anisotropic, via the direction of the hard momentum Q. The 
selection rules conserving angular momentum will come from the integrals. 

3.2 Gluon Exchange Kernel 

Consider the simplest one-gluon approximation to the kernel, H = ATiCpots/ q^'-, 
xP^ — x'P^' + kj< — kj^ ( Fig 1). This is written out as 

H {k — k'rp] X , x' , Q) = AnCpois 



xx'Q^ + (/ct - k'j)'^ 
H{b;x,x',Q) = 87T^CFasKo{Jxx'QW). (5) 



When integrated with non-singular functions of b and x, the dominant region 
is not determined by Eq. ^ but instead the Bessel function restricts to 

Vxx'Qb < 1. (6) 

Clearly Eq. IHlis more accurate than Eq. ^ because the partons entering the 
reaction carry xP, rather than P and can only be scattered through xx'Q"^ 
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momentum transfer-squared. This has been recognized for a long time, and 
common wisdom ascribes an average value of x ~ 1/3 for valence quarks. 
It is commonly accepted that consistent appearance of xx'Q"^ postpones any 
onset of short- distance dominance to higher values of compared to the 
naive implications of Eq. ^ 

A Physics Question: We pause to question how the "uncertainty prin- 
ciple" of EqUl h ~ l/<5, could have misled the field. First, it was not the 
uncertainty principle, but a guess at dominant regions. Somehow x, from the 
longitudinal coordinate, got into a relation between transverse things. We 
seek a physical explanation. 

Suppose we stand near a fast moving charge with energy mass m, at 
impact parameter h. Relativity predicts a pulse of fields with a time scale 
At ~ &/7-^ Since = E/m the fields depend on Atj/b. The partons 
are spread over = xE ~ 1/At. Therefore the fields depend on [mxb). 
This explains how the longitudinal fraction inserted itself. In a perturbative 
transition amplitude between x, x' values, dependence on a/ xx'Qb is generic, 
and nearly kinematic. The time-scale dependence on impact parameter has 
been identified as important. The time-scale (x) dependence is of course 
set by the non-perturbative part of the problem, the x dependence of wave 
functions. 

3.3 Time Scale Smearing: Interplay of x and b 

Exploring the integrals, the 6'™' ~ l/Q'™' counting will still occur, from di- 
mensional analysis, after doing the b integrations in each x, x' integration 
bin. It looks like Eq. ^ and its conclusions will win after all. What is quite 
surprising, and highlighted here, is that there is no simple rule after the 
X, x' integrations are done. The scaling indicated by dimensional analysis 
is "erased" by the x-integrals under broad conditions. The physical origin, 
of course, is that the dynamical time scale distribution in the wave func- 
tions (x distributions) is not something we are priviledged to predict. The 
x-distribution is set by the proton itself in the quiet of vacuum over infinite 
time. The "time-scale smearing" destroys naive use of the uncertainty prin- 

■^This time scale is also a Lorcntz-contracted pancake longitudinal distance scale. Our 
use of "time-scale" is consistent with the impulse approximation. 
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ciple, and its counting of powers. However it does not destroy the use of 
pQCD, which always contains integration over x. 



3.4 The Dominant Power of x 

It suffices to use the kernel Eq. to show the effects. Since our focus 
is the orbital angular momentum, we consider the integrals with factors of 
]ji\m\ gxplicit. The selection rules from the angular integrals e"**^™"*" ""^ ^'^ 
are obvious. We use a Gaussian function <l>(6, x, () = 6^e~''^/'^^"^^0(a:, C) to 
represent the wave functions cutting off large b, as well as a Sudakov model, 
described below. The factor of 6^ is the phase-space to find A quarks close 
together from naive quark-counting. We change variables for the longitudinal 
fractions to 

JU V JiJ JO 1 



We parameterize 



C = x/x' . 



0„(x,c)~x'^+'(i-xrv(c) 



as x — > 1: the C dependence can be left unspecified. For later use parameter 
r is called the "dominant power of x". In discussing the region x ^ 1, we 
can expand (pmi^iC) ^ power series. Knowing the dominant power is all 
that is needed for the arguments to go through. 

3.5 Mellin Method 

We study the large-Q asymptotics by calculating a Mellin transform F{N) 
conjugate to Q^: 



F{N) = r^QMQ) 
Jo (J 



Q 

oo , , . , /•! r roo dQ 



JO JO J JO U 



10 Jo J Jo Q 

The Q integral is carried out easily, 

dQ 



I ^ Q^KoixbQ) = 2-^+^ (xb)-^ [r (^)] ' . (8) 
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Note that emerges just as expected from dimensional analysis. If we 
stopped here, then 6™ would be suppressed by Q~"^. 
The other integrals are then done as 

/" db t dx - a;)'-+ie-^'/(2"') = 

Jo Jo 

^^+A-iv+2 [r(N/2)f r( ^+"^y-^ ) r(2 + r) r(2 - TV + r) 

r(4-iV + 2r) ^ 
We invert the Mellin transform with a contour integral: 

1 pioo 

F{Q) = 7r^ diVF(iV)Q-^. (10) 

Am J-ioo 

The contour in the complex N plane runs in the strip of real-A?^ where F[N) 
converges. Deforming the contour about singularities one-by-one generates 
an asymptotic series in Q. We use notation F{Q) < to indicate depen- 
dence falling at least as fast as Q~^. 



3.6 HHC as Inequality 

Dependence at large Q now comes from singularities to the right of the con- 
vergence strip. There are two distinct types. 

The singularities of r( 2+^+^--^ ) are simple poles at 

N =2 + m + A + 2K, K = 0,1,2... 

These are exactly the singularities creating the naive power counting of 
HHC. The existence of these singularities implies 

F{Q) KQ-"^-^-^ {hhc region only), (11) 

from these singularities alone, and barring a zero factor canceling the poles. 
As expected 6"* here the HHC results are reproduced. 

The dominant power r locates the other singularities, seen in the pole 
from r(2 — N + r). These poles are independent of the 6™ dependence. 
This implies an additional contribution to power-behavior, not regulated by 
bQ < 1. The poles and power behavior are 

N = 2 + r + K K = 0,1,2... 
F{Q) ^ Q~'^~^ , {dominant power region) (12) 
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regardless of the value of m. If such poles dominate, there is no power 
suppression of 0AM m 7^ 0. 

As a resuh we have the following inequality for the large Q dependence 
of the form factor, 

nQ) < ^ 

c = Min{r, m + A). (13) 
There is no power suppression of 0AM m 7^ in general. 

Rule of Thumb: We have extracted an improved HHC-rxde of thumb: the 
asymptotic dependence is determined by the minimum of powers m, r in 
the integrands. The rule has been deliberately simplified to powers, suppress- 
ing logarithmic factors readily calculated, by further evaluating the Mellin 
inversion. These logarithms are obtained from the Mellin series and are 
a separate phenomenon from the asymptotic-freedom logarithms of pQCD. 
The origin of the inequality Eq. is time-scale smearing: due to it, the 
suppression of 0AM expected by naive counting is transferred to smearing 
inside the longitudinal wave functions, and helicity-flips can readily compete 
with helicity non-flip. 

Effects of the Large b Regions The singularities to the left of the con- 
vergence strip determine a series for small Q ^ 0. These singularities are 
traced to the Gaussian model, Sudakov 6-cutoff, or the large b cutoff. We are 
not concerned with these singularities, showing that the large b cutoff drops 
out of the power-laws stated above. However the numerical normalization 
of integrals depends on the cutoff method. It is a separate question needing 
separate analysis to find the numerical fit of models to data. Our numerical 
work is presented momentarily. 

On this reasoning there is always a part of the calculation which is 
"strictly perturbative" . We may calculate the large b regions only with lim- 
ited reliability: fine! The small b regions known to be calculable behave as 
we claim, and are the only regions under conceptual dispute. 

3.7 Discussion 

As far as we know this is the first time that HHC has been observed to fail for 
large form factors. Yet previous work has made related observations that 



11 



the powers of x can determine the dependence. Feynman's mechanism 
[77j concentrated entirely on the region x — 1 and ignored b. Feynman's 
work predated pQCD, and there are certainly regions in the b integrals of 
that mechanism that would not be calculable, so the model is moot. Soper's 
1977 paper has a clear statement that the minimum power of x or kr 
rules the results. The context there was the Drell-Yan-West [2H1 relation, 
which connects the x —>■ 1 behavior of inelastic structure functions with the 

dependence of the form factor. Soper shows that one cannot prove the 
Drell-Yan-West relation deductively, but one can get an inequality and force 
the relation by choosing wave functions. 

The literature is clouded here by treatments assuming that the two regions 
of X ^ 1 and 6 — must give the same scaling, and forcing a result by 
circular logic. As we have just shown, the regions are different and no general 
rule can be made. Underlying this is the fact of perturbative calculability 
of short-transverse distance not being on the same footing as perturbative 
models of the x —>■ 1 dependence. For instance Brodsky and Lepage ^] 
discuss the endpoint x —>■ 1 contribution of the proton form factor Fi. The 
authors argue that for the limit Xi —>■ 1, pQCD implies an m = wave 
function of order a'j. perturbatively calculated to go like (1 — xi)^. On this 
basis they find that a contribution to Fi scaling like 1/Q^ independent of 
the powers of b inside the integrals. The Q dependence of this contribution 
is same as that obtained from the b ^ 1/Q region. Although perturbative 
analysis is not a valid approach to calculating wave functions, the result is 
consistent with ours. 

Our results are due to time-scale smearing and should not be attributed 
to "end-point singularities" as the term is commonly used. Inside the hard 
scattering kernels in pQCD are combinations of inverse powers of x or 1 — x. 
We are not exploiting these inverse powers in the evaluation of the Mellin 
moments. Our basic point is that the x — > 1 limit of the m ^ wave 
function is not known. This uncertainty leads to the inequality given in Eq. 
IT^ The end point singularities, when present, further enhance the end point 
contributions and expand the range of allowed values of r, the dominant 
power of X, for which there is no power suppression of 0AM. We are also not 
concerned with endpoint wave functions (a la CZ) for the same reason: their 
relevance or lack of it is another subject. 

To summarize the logic so far, for broad classes of x dependence of the 
wave functions, which is unknown, the non-rule of HHC is transformed to a 
rule of HHNC, hadron helicity non-conservation. This is a new asymptotic 
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prediction, and proves that Q F2(Q^)/Fi(Q^) ~ const can be an outcome of 
the theory up to the highest Q. To be equally fair, different broad classes 
of wave functions give Q'^ F2{Q'^) / Fi{Q'^) ~ const. Our primary accomplish- 
ment enlarges the sphere of allowed possibilities. Measurements are still 
needed to determine what protons are. 



4 Numerical Studies, Sub- Asymptotic 

The scaling rule Q F2{Q'^) / FiiQ"^) ~ const describes experiments at labora- 
tory Q-values very far from asymptotic. We explored this region numerically. 



4.1 Studies with the Pion 

The pion provides a simple test system. We test the dominant integration 
regions by inserting factors of 6"^ in the integrands, corresponding to m-units 
of 0AM. The moment < h{Q) >t^ is defined by 

, / dxdx' cPhhJ^JQ , x, x', h) ^. 

"^'^^ w) ■ 

The pion form factor Ft^{Q'^) is given by 

F^(q2) ^ J dxdx'd%J^^{Q, X, x\ h) (15) 

with the kernel Tt,{Q, x, x', h) defined by 

J^^(Q, X, x', b) = j dxdx'dbb(j){x')as{iJ,)e-^^''''''''''^'>Ko{V^bQ)(f){x)<l>{b) . 

(16) 

Here S{x, x', b, Q) is the Sudakov form factor. A model soft wavefunction 
$(6) = exp(— 6^/2a^) is also included, where the parameter a = 1/Aqcd- 

In Fig. ^ we plot the moment < b{Q) using a wave function of the 
form 

0(x) ~ [X{l-X)]\ 

We study values of 6 = 0.2, 1, with and without the Sudakov form factor. 
In the latter case the momentum scale /i = Q/A of the strong coupling is 
imposed, and the b integrals are cut off at 6 = 1/Aqcd- From the Figure it 
is clear that for < 100 GeV^ the moment falls much slower than the asd 
prediction of l/Q, regardless of the x-dependence of the wave function used. 
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Figure 1: The moment < b{q) >t, of the pion form factor kernel, as defined 
in the text, using the [x{\ — x)Y form of wave function for 5 = 0.2, 1. The 
moment decreases with Q much slower than the 1/Q behaviour expected 
in the asd HHC model. Results are shown with and without including the 
Sudakov effects. The sohd (l/Q^'^^) and the dashed (l/Q^'^) hne represent 
a simple power law fit at small Q'^ for 5 = 1 and S — 0.2 respectively. 
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Figure 2: The proton form factor Fi for moderate Q^, extracted using the 
Jlab data for Ge/Gm and the SLAC data for Gm- The ratio Ge/Gm is 
obtained from the parameterization P] jJ^p^^ = 1 — 0.13((5^ — 0.04). The 
sohd hne represents Fi ~ 

4.2 Studies with the Proton 

We turn to the proton form factors. We note that the JLAB data for QF^jF^ 
is flat even below the Q-range where Fx ~ begins to fit. In Fig. |21we 
show a plot of Q^Fi. The solid line in the plot corresponds to the behaviour 
Fi ~ It is clear from the plot that the scaling Fi ~ l/Q^ is seen only 

for > 6 GeV^, which is larger than the momentum regime explored at 
JLAB so far. This is cause for concern. It is nevertheless entirely possible 
that the scaling observed in the ratio is not overly sensitive to Q^, and will 
continue to larger Q^. We investigate this in greater detail numerically. 

To probe the dominant integration regions, we turn to calculating b- mo- 
ments of the proton form factor kernel. These are multidimensional inte- 
grals of which our analysis in the previous section determines but a low- 
dimensional strip. We spare the reader a listing of dozens of Feynman dia- 
grams and substantially complicated kernels listed in the literature The 
form factor Fi in impact parameter coordinates can be written symbolically 
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Figure 3: The moment of the transverse separation bi and 62 for proton 
using the wave function of the form {xiX2X3)'^-^ . The < 62 >p moment is 
shown both with and without including the sudakov form factor. The small 
dashed (l/Q^-^^), medium dashed (l/Q^-^) and the large dashed 
lines represent simple power law fits at large Q^. 
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as [2S| 



~[ 27 Jo Jo Jo 

xHj{xi, x[, bi, Q, tji,tj2) "^jixi, x'i, w) 
xexp[-S{xi,Xi,w,Q,tji,tj2)] , (17) 

with (dx) = dxidx2dx36{l — Xi — X2 — x^). The variable 9 is the angle 
between bi and b2 and Xi and x'^ refer to the initial and final x variables. 
The expressions for the hard scattering Hj, the Sudakov form factor S, and 
function are given in Ref. Now defining 



Fi{Q) = J hidhh2dh2Tp (18) 
from Eq. [TTl we define moments < hj{Q) >p as follows: 

<w)>P p^^Q2) ■ uyj 

The function is where the linear combinations of the products of initial 
and final x wave functions are found. The most singular part of the kernel 
in the limit a;i — > 1 and x'l — > 1 is obtained from the i^i^E'i term, which is of 
the form, 

Ko U{1 - xi)il - x\)Qh^ Ko (T^gfes) 0(x.)0(x9 

Hl^l ~ ; 7- (20) 

(l-xi)(l-x'i) ^ ' 

Here Kq is the modified Bessel function of order zero. 
For the test we explore a wave function (j){xi) given by 

(f){Xi) ~ {XiX2XsY . 

The numerator in is then proportional to ( ) . The Bessel 

functions imply that in the limit of large Q, (1 — Xi)(l — x[)b1 l/Q"^ and 
X2X2b2 —>■ 1/Q^- As long as 5 < 0.5 the time-scale smearing will dominate. 
For 6 = 0.5 we get Fi 1/Q^ even though a dominant region in bi and 62 
is independent of Q ! For 6 < 0.5 the moments of 61 and 62 should also have 
a region independent of Q. 



17 



We check these predictions by performing the calculation using S = 0.5. 
We ignore the Sudakov form factor for this test, because it is a side issue. 
We evaluate the strong coupling at In figure El we plot the moment 

of the transverse separations bi{Q) and 62 (Q) as a function of Q. We find 
that for Q"^ > 100, the moments < 62 >p and < bi >p fall asymptotically as 
l/Q^'"^ and respectively, in agreement with our analytic expectations. 

If we include the Sudakov form factor the moment < 62 >p falls as 
asymptotically. This is only a slightly stronger decay compared to the earlier 
case. In contrast to earlier expectations, the Sudakov form factor does not 
much suppress the importance of the end-point region. 

The results with end point dominated COZ [221 ED] x-dependence are 
similar. In this shown in Fig. El we find that the < 62 >p decays 

very slowly with Q for a wide range of Q. This is quite interesting, because, 
as shown in Fig. the form factor Fi itself does not show good scaling 
in this momentum regime, if the sudakov form factor is not included. The 
calculation indicates that the scaling seen in the moment, and hence the ratio 
QF2/F1, is more general than that seen in Fi. This is not totally unexpected 
in a ratio. At very large Q it starts to fall faster, but only as well 
below the supposed 1/Q rule of asd assumptions. 

4.3 The Ratio Q F2{Q^) / Fi{Q^) 

Based on these studies we come to a prediction for the ratio (5-^2 (<5^) / -^i(Q^)- 
This prediction is based on power-counting, and the dominant integration re- 
gions probed by the moments. It is as complete a prediction as now possible, 
taking into account that one unit Am = 1 of orbital angular momentum is 
needed for F2 to proceed for » GeV^. Up to small and model dependent 
corrections, of order Q^-^^ , the power counting gives the scaling behavior but 
not the normalization of the data, via the relation 

< b2iQ) >p~ ^/'^^P ~ constant. (21) 

What if this prediction fails? It can fail, according to our asymptotic 
studies, if the dominant power r is "large" . Large studies of F2 not only 
probe quark 0AM, but also they can tell us details of the dominat x power 
associated with m 7^ 0. 
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Figure 4: The moment of the transverse separation < 62 (Q) >p for the proton 
form factor kernel using the COZ wave function [SH]- The sohd and 
dashed hues represent simple power law fits at small and large 

respectively. The dependence on Q is much weaker than predicted by asd 
relations, and supports a flat prediction for QF2/F1 at JLAB momentum 
transfers 
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Figure 5: The dependence of the calculated proton form factor Fi us- 
ing the COZ wave function [211 • Results are shown both with and without 
including the Sudakov form factor. Details are in the text. 



4.4 Positivity of Ge 

In the Appendix of a 1959 Reviews of Modern Physics article there are 
listed two form factors denoted A and B, now called Ge and Gm, which are 
linear combinations of Rosenbluth's: 

Gm = Fi + KF2] 
Ge = F,- KTF2; 



4m^ 



(22) 



and K is the anomalous magnetic moment. For some reason these are called 
the Sachs form factors [S2]- Yennie et al gave |3T] the alternative linear 
combinations simply to emphasize that the definitions of form factors was 
arbitrary. 

It is interesting to observe that if the trend of QF2{Q'^) / Fi{Q'^) ~ const 
continues, the value of G^; will reach zero, and cross to a negative value. 
One wonders whether there is a physical significance and a barrier to 
this unusual occurence. The meaning of G^; is angular momentum AJ^ = 
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spin-preservation in the Briet frame. If Ge = the proton spin must flip 
in the scattering, which seems industrially useful. We find nothing special 
about this. The amplitudes in the spacelike region cannot be limited further 
than general principles of analyticity, Lorentz and gauge invariance, and so 
on. It is perfectly consistent to arrange timelike discontinuities so that the 
sign change of occurs without violating anything holy. 

5 Concluding Remarks 

Our calculations show that the (newly revised) asymptotic behaviour of pow- 
ers of b is achieved only at very large Q^. Physics has many asymptotic 
predictions, which by construction have suppressed information needed to 
know when they will apply. To repeat: methods inherent in asymptotic pre- 
diction strongly tend not to tell you where the prediction will be valid. In 
QCD it was early thought the "asymptotic" regime had to occur at ener- 
gies very large compared to Aqcd, the strong scale. This was built wrongly 
into dogma. Very often asymptotically large logarithms are needed to justify 
some approximation. We have shown here that the scales of asymptopia for 
form factors are vastly beyond experimental comparisons. Asymptopia is 
effectively meaningless, because all the new physics of very high energies has 
been left out. 

Meanwhile, the finite effects of 0AM are hardly suppressed at all. 
They are suppressed in explicit calculations by even less than the revised 
asymptotic behavior. The pyramid of assumptions that the s-wave distribu- 
tion amplitudes are meaningful falls into very grave doubt. To put this more 
directly, we don't have a reason to use a distribution amplitude any more. 

The usual approach to pQCD, in which wave functions are assumed to 
be unknown, predicts HHNC: hadron-helicity non-conservation. There is an 
inequality that helicity-fiip processes do not lead but can have equal power 
with helicity non-flip processes. 

The experimental observation of R{Q) = QF2/F1 appears to have a 
dual meaning. The ratio is a very robust quantity, which remains flat even 
in the regime where Fi is not clearly dominated by quark-counting. The 
same ratio is an important asymptotic quantity, which says things about the 
wave functions. The common notion that proton wave functions are "cubic" , 

'^Quantum computing, which often invokes high-speed, high-precision spin-flips, may 
or may not ever need the physics of multi-GeV spin flips, but it stands ready. 
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namely going like XiX2Xs, is ruled out for wave functions calculating F2. 
Applications of the CZ wave functions are ruled out for F2, if the ratio R 
continues to be fiat at large Q^. 

Predictions for protons are mirrored in predictions of the same type for 
neutrons. If the flat R{Q) ratio is followed for neutrons, then the hallowed 
Galster fits to neutron data will eventually fail [23] ■ We predict a fiat R 
ratio for neutrons on the basis of isospin independence of QCD. 

The study of polarization transfer in large nuclei, e{A, e')j? should yield 
further information. If the existing form factors are dominated by asd, which 
we do not believe but is still worth testing, then nuclear filtering will not 
make the distances shorter, and the ratio R should be fiat with A >> 1. 
The regime of small A ~ 10 has no advantages and many complications due 
to few-body effects, and predictions are more difficult. If the existing form 
factors are dominated by quark mass effects, then nuclear filtering will not 
make any difference and the ratio R should again be flat. If quark 0AM 
is responsible for the flatness of R as we believe, then flltering will kill the 
large transverse extent of large 0AM jSHl ISH] , F2 should be depleted relative 
to Fi, and R will decrease with A >> 1. We intend to dedicate a study to 
quantifying these predictions. 

JLAB has made a pivotal experimental discovery which will be a perma- 
nent subject of discussion. The experimentally observed flatness of QF2/F1 
is a signal of substantial quark orbital angular momentum in the proton. 
Higher momentum transfer measurements would be helpful in conflrming 
this interpretation. The numerical size of R cannot be converted directly to 
a wave function, because it is only a single number. But the value of R can 
rule out models which omit quark 0AM. Further studies at higher may 
separate constituent quark models with 0AM jB] , which tend to have a scale 
(the quark mass) forcing turn-over of R with increasing . 
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